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We present a theoretical approach to determine the electronic properties of nanoscale systems ex- 
hibiting strong electron-electron and electron-phonon interactions and coupled to metallic electrodes. 
This approach is based on an interpolative ansatz for the electronic self-energy which becomes exact 
both in the limit of weak and strong coupling to the electrodes. The method provides a generalization 
of previous interpolative schemes which have been applied to the purely electronic case extensively. 
As a test case we consider the single level Anderson-Holstein model. The results obtained with 
the interpolative ansatz are in good agreement with existing data from Numerical Renormalization 
Group calculations. We also check our results by considering the case of the electrodes represented 
by a few discrete levels which can be diagonalized exactly. The approximation describes properly 
the transition from the Kondo regime where electron-electron interactions dominate to the polaronic 
case characterized by a strong electron-phonon interaction. 
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I. INTRODUCTION 

While Kondo physics due to electron-electron in- 
teractions have received a long standing attentionii^, 
the electron phonon interaction and its interplay with 
electron-electron correlation effects have only attracted 
more recently the interest of researchers in the field 
of nanoscience^ . This is probably due to the impor- 
tant role played by both electron-electron and electron 
phonon interactions in the transport properties of small 
molecules^i^i^. In this respect, the Anderson-Holstein 
model (AH)ii2ii2ii2iiii, including these correlations in a 
local way, is a paradigmatic simple model to study the 
physics of these systems. 

Several theoretical methods have been success- 
fully developed for describing Kondo correlations in 
the single level Anderson model like the Numerical 
Renormalization Group (NRG) ^^d^ , Quantum Monte 
Carlo (QMC)ii, the Bethe Ansatai^ or perturbative 
methodsi^iii. However, the extension of these methods 
to deal with both electron-electron and electron-phonon 
interactions has certain limitations. We mention here 
the NRQi^iiai^OiSi and QMC methods^ that have been 
recently applied to the Anderson-Holstein Hamiltonian: 
in the NRG approach it is difficult to reach the strong 
polaronic regime in which the electronic density spec- 
trum develops a very rich structure of phonon induced 
peaks over a broad range of energy; in the QMC ap- 
proach, it is also difficult to resolve these features due 
to the limitation of using a minimum temperature in the 
calculations. On the other hand it would be extremely 
demanding from a computational point of view to imple- 
ment these numerically exact methods in more complex 
models describing realistically systems like molecules or 
Quantum Dots (QD) coupled to electrodes. 

These limitations suggest the convenience of looking 
for other alternative methods that would yield approx- 



imate but sufficiently accurate solutions for these kind 
of systems for a broad range of parameters. In the 
case of purely electronic models an approach that has 
proved to be very useful is the interpolative self-energy 
method^iSiiS^iS^iS. This approach is based in the anal- 
ysis of the properties of the exact self-energy in the weak 
(perturbative) and strong coupling limits which allows to 
define an approximate and simple self-energy interpolat- 
ing appropriately between the two limits. This scheme, 
originally derived for the equilibrium Anderson model, 
has been generalized and extensively used in different 
contexts: multilevel QDsS^, out of equilibrium transport 
through a single leve P^'^° and in combination with Dy- 
namical Mean Field Theory (DMFT)^i2£, to analyze the 
Mott transition in Hubbard- like models. Its potential to 
be extended to more complex situations makes this ap- 
proach particularly useful. The aim of this work is to 
present a generalization of this interpolative approach to 
deal with both electron-phonon and electron-electron in- 
teractions simultaneously. 



In order to introduce the method we will consider in 
this paper the Anderson-Holstein model which includes 
both type of interactions. We first discuss the spinless 
case which is instructive as will allow us to illustrate 
how the interpolative arguments also work for the purely 
electron-phonon interaction case. This case will be ana- 
lyzed in section II. The extension to the case where both 
interactions are present will be addressed in section HI. 
In both cases the accuracy of the method is tested by 
comparing the results for the spectral density with ex- 
act numerical results for finite systems and with existing 
NRG calculations in the literature. Finally, in section IV 
we present a discussion and the conclusions of our work. 
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II. INTERPOLATIVE SELF-ENERGY FOR THE 
SPINLESS ANDERSON-HOLSTEIN MODEL 

We consider first the simple spinless Anderson-Holstein 
Hamiltonian describing a single non-degenerate elec- 
tronic level, eoj coupled linearly to a local phonon mode 
of frequency uja and to an electronic reservoir represented 
by a continuum of states of energy e^. 

k 

+ T.{^kic + c.c) (1) 

k 

where Vj^ is the matrix element between the reservoir 
states and the localized level, fi is the level occupation 
operator and A the electron-phonon coupling parameter. 
It should be noticed that in spite of being simpler than 
the spin dependent case, this model has not an exact so- 
lution for an arbitrary level occupancy. Only in the limit 
of a completely empty (or full) level, an exact solution is 
availableS-. 

As in the purely electronic case^^i^ the interpolative 
approach will be based on a special property of the elec- 
tronic self-energy which exhibits the same mathematical 
form when expanded in the interaction parameter both 
in the atomic (Vk — > 0) and in the perturbative limit 
(which in this case corresponds to A — > 0). 

We first consider the atomic limit, — > 0. In this 
case Eq. ([T]) can be diagonalized by means of a canon- 
ical transformation (see Refs,'^^''^^). The electron Green 
function can then be straightfordwardly obtained: 
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where io — to ~ A^/cjq and (n) is the level occupation. 

This expression can also be written as a continuous 
fraction which can be useful for computational purposes 
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The electronic self-energy in the atomic limit can now 
be obtained from the Dyson equation S^"*) — lo ~ en ~ 
Qiat)-i -y^rhere en = to — 2(A^/ijJo)('^) is the energy level 
corrected by a constant term which can be identified as 
the Hartree contribution. In the limit of small coupling 
X/uJo << Ij and up to order A^, this expression tends to 



E('^*)(a;) « A^ 



1- < n > 



< n> 



(4) 
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FIG. 1: Lower order diagrams contributing to the proper self- 
energy. Diagrams a) and b) correspond to the spinless case. 
Diagrams c), d), e) and f) to the spin dependent Anderson- 
Holstein model. 



On the other hand, the self-energy of the model can be 
calculated also up to order A^ by means of perturbation 
theory in the electron-phonon interaction. The two dia- 
grams contributing to order A^ are depicted in Fig. la) 
and Fig. lb). Diagram a) gives simply the constant 
Hartree contribution —2{\^/tjjQ){n)f^, where (n)^ is the 
unperturbed level occupation for an effective level , ee//, 
which will be fixed by an appropriate self-consistency 
condition as commented below. 

The correlated part of the self-energy to order A^ is 
then given by the diagram of Fig. lb and has the expres- 
sion: 



de- 



p(0)(6) 



de- 



3(0) (e) 



LO — e + ujQ + if] 



(5) 



where p^°\lo) = F/ [(w - ee//)^ -I- F^] /tt is the level 
density of states of the one-electron unperturbed case, 
F = TrJ2k \^k\'^S{LO — tk) which is taken as constant and 
/i is the reservoir chemical potential. 

In the limit F — ^ the above expression tends to 



1- < n >o 



<n>o 



^F{lo) 



- Eeff + UJQ ^ 

(6) 

By comparing Eqs. (H]) and ^ one can notice that 
the exact atomic self-energy up to order A^ and the sec- 
ond order self-energy in the limit F — > have the same 
functional dependence on the occupations and frequency, 
as was the case in the purely electronic case. Therefore 
it is possible to find an interpolation following the lines 
discussed in Refs.— the only difference being the more 
complex expression for the atomic self-energy due to the 
multi-phonon satellite structure. The interpolative self- 
energy is defined by the following ansatz: 
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(7) 



where is the inverse function defined by Eq. ([6]). 
It is straightforward to check that this ansatz interpo- 
lates correctly between the weak and strong coupling 



limits. The atomic limit is recovered because in the 
limit r/A ^ 0, ^ ui. On the other 

hand, in the limit of small electron-phonon coupling, 
j.(at) 1^-1 [s(2){t^)]} ^ S(2)(tj), and the results of per- 
turbation theory are recovered. ^From Eq. ([6]), we 
can straightfordwardly calculate F^^ [S*^^^(a;)] = Q{lj) 
which is given by the roots of a quadratic equation: 



2E(2) 



■(1 - 2 < n >o) + 



A2 



(8) 



where the -I- sign in front of the square root has been 
chosen to yield il(a; — > ±oo) = uj; this guarantees the 
recovery of the atomic limit for if F/A ^ 0. More- 

over, we have checked that Eq. ([5]) defines an analytical 
function of lo with poles in the positive complex half- 
plane: this is a necessary condition to have a self-energy 
with the correct analytical properties. 

We now discuss how to determine Ce// by imposing 
an appropriate self-consistency condition. Several pro- 
cedures have been proposed in the literature: The sim- 
plest one is to impose charge consistencj*2^i2ii2S between 
the one-electron and interacting problems, (n) — (ti)q. 
As discussed in^^ this procedure only fulfills the Friedel 
sum-rule in an approximate way although the agreement 
is rather good for a broad range of parameters. An alter- 
native possibility is to fix the effective one-electron level 
by imposing the fulfillment of the Friedel sum rule^^ . The 
self-consistency condition is then given by the equation: 

<n>=- arctan ( p ) 

where the occupation (n) is calculated from the interact- 
ing retarded propagator G^{uj) = [lj — €h ~ E(aj)]~-^: 

<n>= / dujlmCiuj) (10) 

^ J~oo 

This last condition seems to give a better description of 
the density of states around the Fermi level, specially in 
cases with a large electron-hole asymmetry. 

This interpolative ansatz reproduces very accurately 
both the limit of weak and strong coupling regimes while 
giving a good quantitative description in the intermedi- 
ate range. In order to check this approximation we first 
compare its predictions for a discrete version of the spin- 
less Anderson-Holstein Hamiltonian in which the elec- 
trode is replaced by a three site tight-binding chain which 
can be numerically diagonalized. This is a useful model 
which allows a detailed comparison with exact results be- 
tween the weak and strong coupling regimes. For a given 
value of the electron-phonon coupling A, the configura- 
tion Hilbert space is truncated in the number of phonon 
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FIG. 2: (Color online) Localized level spectral density for a 
discrete spinless Anderson-Holstein model with an electrode 
consisting of a three atom tight-binding chain for different val- 
ues of the electron-phonon coupling parameter A with ujo = 1. 
From upper to lower panel, A = 0.6, 1, 1.4. The (blue) contin- 
uous line corresponds to the interpolative solution while the 
(black) dotted one to the numerical diagonalization results. 
The diagonal energy levels in the chain are taken as zero and 
the first-neighbors hopping parameter is fixed at t = 1. The 
hopping parameter between the dot level and the chain is 
V = 0.5 and the level position is the one corresponding to an 
e-h symmetric case. A small imaginary part, r) = 0.05, is in- 
cluded in both spectral densities to facilitate the comparison. 
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modes in such a way that convergence in the level spec- 
tral density is achieved. 

In Fig. [2] we compare the level spectral density, 
Pq{uj) = — ImG"'(w)/7r, obtained from the numerical di- 
agonalization with the approximate solution given by the 
interpolative self-energy for a case with electron-hole (e- 
h) symmetry. As shown in the figure, the interpolative 
approximation yields in all the cases an excellent approx- 
imation to the exact DOS. Notice that the range of pa- 
rameters goes from the weak coupling limit (A = 0.6), to 
the intermediate (A — 1) and the strong coupling case 
(A = 1.4). For lower and higher values of A both spectra 
are practically indistinguishable. 

Results for the continuous model corresponding to 
Hamiltonian Eq. ([T]) are shown in Fig. [3l As in Fig. [2] 
the three cases go from the weak to the strong coupling 
limit showing the interpolative solution and the results 
provided by second-order perturbation theory in the self- 
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FIG. 3: (Color online) Localized level spectral density for 
the continuous spinless Anderson-Holstein model for differ- 
ent values of the electron-phonon coupling parameter for an 
e-h symmetric case with F — 0.25 and ujo — 1 (from upper 
to lower panel, A = 0.3, 1, 1.4). The (blue) continuous line 
corresponds to the interpolative solution and the (black) dot- 
ted line in the upper panel to the second-order self-energy 
approximation. 
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FIG. 4: (Color online) Localized level spectral density for a 
discrete spinless Anderson-Holstein model with an electrode 
consisting of a three atom tight-binding chain for an asymmet- 
ric electron-hole case. The (blue) continuous line corresponds 
to the interpolative solution while the (black) dotted one to 
the numerical diagonalization results. The renormalized level 
position is eo = —0.2 with A = 1 and wo = 1- The diagonal 
energy levels in the chain are taken as zero and the first- 
neighbors hopping parameter in the chain is fixed at t = 1. 
The hopping parameter between the dot level and the chain 
is V = 0.25. A small imaginary part, r] = 0.05, is included in 
both spectral densities to facilitate the comparison. 



energy. It is interesting to realize that our interpolative 
solution is a fair approximation of the second order per- 
turbation theory for small A; however, for the interme- 
diate or the strong coupling limit, our solution develops 
the typical structure of the multiphonon excitation spec- 
trum which is more clearly seen at high A. In addition to 
this multi-phonon features, the DOS exhibits a narrow- 
ing resonance at the Fermi level when approaching the 
strong coupling limit. It is straightforward to calculate 
the weight of this resonance from the value of the self- 
energy at the Fermi level, giving an exponential decaying 
law « exp{— X'^ / ujq) . 

For the sake of completeness we show in Figs. |4]and[5] 
results for the level DOS for an electron-hole asymmetric 
case. In Fig. 2] we compare our interpolative solution 
with exact diagonalization results for the same discrete 
model used in Fig. [51 Notice the excellent agreement 
between both solutions which correspond to a level oc- 
cupancy (ridiag) ~ {n) ~ 0.83. In Fig. [5] we show the 
local DOS for the full continuous asymmetric model for 
a case with a level occupancy (n) = 0.82 (see caption 
for details). Notice how the non-interacting resonance 
(blue dashed line) having a Lorentzian-like shape around 
eo is changed into a density of states showing a satellite 
structure associated with phonons. 

As a final remark, we would like to point out that the 
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FIG. 5: (Color online) Localized level spectral density for the 
continuous spinless Anderson-Holstein model for an electron- 
hole asymmetric case. The parameters values are A = 1, 
r = 0.25, Luo ~ i and eo = —0.25. The (blue) continuous 
line corresponds to the interpolative solution and the (black) 
dashed line to the non-interacting case with eo = eo- 

interpolative ansatz recovers the exact solution in the 
limit of a fully empty or fully occupied level. In this 



case Q ^ Lo + iT and the exact propagator in this limit, 
G{lu) = G^''*) {ui + iT) is recovered^.. 



III. THE SPIN DEGENERATE 
ANDERSON-HOLSTEIN HAMILTONIAN 

In this section we will show how this interpolative 
scheme can be extended to the presence of both electron- 
electron and electron-phonon interactions, by considering 
the spin degenerate Anderson-Holstein Hamiltonian: 

+ A(5t -f 6) ^n„ + ^(^^cj^c. -fee.) (11) 

ka- 

where U is the Coulomb interaction in the localized level. 
The atomic Green function can also be obtained from 
a canonical transformation of this Hamiltonian in the 
limit Vfc — !■ giving a renormalization of both the level 
position Iq = eQ — X^/ujq, and the Coulomb interaction 
U = U — 2\'^ /ujo- As a function of these parameters the 
atomic Green function takes the formi^: 



m=0 



Lu — eo — mujQ 1^ _ _ (7 _ 

< ncr(l - ng) > ^ < ngng- > 
u! — io + muiQ uj — io — U + mujo 



(12) 



For computational purposes, this expression can also be 
written in a continuous fraction form. Following an anal- 
ogous procedure to the one used in the spinless case, 
the atomic self-energy can be calculated from the above 
Green function as S*^"*^ = uj — tHa — G^'''^~^ where the 
Hartree level has in this case the expression 



^Hcr — Co 



(13) 



Expanding now the correlated part of the atomic self- 
energy up to order and we obtain the equation: 



1- < n<j > 



U 



LO — ~ bjQ w — eo + cjQ 

2 < fig > (1- < fig- >) 



(14) 



On the other hand, the self-energy calculated by means 
of the lowest order diagrams (see Fig. Id) and f))has 



the following expression in the F — > limit (up to the 
Hartree constant contribution given by diagrams 2c) and 
2e)): 



1- < fla- >o 



< ha >0 (1- 



< rig >o 
- ee// + t^o 



<"^=^-^^F.(c.) (15) 



+ 

Notice that when both electron-electron and electron- 
phonon interactions are present, both equations have still 
the same functional dependence. Thus one can introduce 
an interpolative self-energy using the same approach de- 
veloped for the spinless case; in particular we define the 
self-energy by the ansatz: 



(16) 



^ being the inverse function defined by Eq. (|15p . 
These equations define the interpolative self-energy we 
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are looking for, yielding the appropriate behavior for 
small values of A and U, and in the atomic limit, when 
Vk 0. An important difference with the spinless case 
appears because the atomic selfenergy for the Anderson- 
Holstein Hamiltonian depends on the correlation function 
{ncrfia) which has to be calculated in a self-consistent 
way. This is discussed in the appendix. 



In order to obtain F„ 



= flcr{oj) from Eq. 



(fT5)l it is now necessary to solve the following cubic equa- 
tion: 



+ <ng >o (1- < >o)i^l - ^o) 



(17) 



where D,^ ~ — Seff- This equation has two clear 
limiting cases. In the limit 17 >> and U >> X, the 
solution of Eq. pT|) tends to 



, C/2 < 71s >0 (1- < n<T >o) 
^eff H (T^ 



(18) 



which corresponds to the purely electronic case in the 
absence of electron-phonon interaction (see Refsi^i2ii2^). 
On the other hand, in the opposite limit, C/ << and 
U « X, Eq. lfT7|) reduces to a quadratic equation like 
the one in the spinless case and with a similar solution 
(Eq. ([8])) which would correspond to the strong polaronic 
limit: 



)I]i2) ^X^\n„ + (l-2< n„ >o) cool (19) 



Solving the full cubic equation requires choosing the 
root that has the correct physical behavior. As in the 
spinless case this can be achieved by starting from the 
solution that behaves as ^l{uj — > ±oo) k, lo and impos- 
ing continuity for all frequencies. This solution exhibits 
a transition from a regime with [/ > in which 
qualitatively behaves like the one in the electronic case 
to a regime for ?7 < in which ^a{^) rather behaves like 
the one corresponding to the quadratic equation of the 
polaronic limit. The transition takes place around J7 w 0, 
its precise value as a function olX,U and F can be calcu- 
lated by a careful analysis of the analytical properties of 
the roots of Eq. p?)) . We have furthermore checked that 
in the polaronic regime [U « 0), the solutions of Eq. 
(fT7)) are well described by the roots of the much more 
simple quadratic equation: 

{^ll - co2)j.(^2) =\^\n„ + (i-2< >o) c^o 

+ C/^ < ns >o (1- < ns >o)n„ (20) 

which is obtained by neglecting in Eq. p?]) the term 
—U'^ < rig >o (1— < ng >o)uJq. It is easy to check that 
both limits — > and X/U ^ oo are recovered by this 
simpler equation. 




FIG. 6: (Color online) Localized level spectral density for a 
discrete Anderson-Holstein model with the electrode simu- 
lated by a single discrete level (taken as the zero of energy) 
for different values of the electron-phonon coupling A with 
ijo = 1. From upper to lower panel, A = 0.9, 2.6. The (blue) 
continuous line corresponds to the interpolative solution while 
the (black) dotted one to the numerical diagonalization re- 
sults. The hopping parameter between the dot level and the 
chain is V = 0.25 and the dot level position is the one corre- 
sponding to an e-h symmetric case. A small imaginary part, 
r] — 0.1, is included in both spectral densities to facilitate the 
comparison. The Coulomb interaction parameter is U — 8. 



In order to test the interpolative solution presented in 
this section we first compare its predictions with the re- 
sults of exact numerical diagonalizations for a finite sys- 
tem. As the Hilbert space in the case were both electron- 
electron and electron-phonon interactions are present is 
much larger than in the spinless case we take a minimum 
cluster in which the electrode is replaced by a single dis- 
crete level. In Fig. [S] we represent the evolution of the 
level DOS for increasing values of the electron-phonon 
coupling parameter A. It is interesting to notice that 
in spite of the simplicity of this discrete model the local 
DOS exhibits as a function of A a transition from a regime 
where Coulomb correlations dominate (Kondo regime, 
see upper panel) to a strong polaronic regime where the 
effective electron-electron interaction becomes negative 
and is characterized by a gap in the quasi-particle spec- 
trum and the appearance of well defined multi-phonon 
satellites (see lower panel). As can be observed the inter- 
polative solution is very close to the exact diagonalization 
results throughout these regimes. 

In Fig. [7] we show our calculations for the local- 
ized level spectral density for the full Anderson-Holstein 
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IV. CONCLUSIONS 



20- 




FIG. 7: Localized level spectral density of the continuous An- 
derson Holstein model for a fixed value of the Coulomb in- 
teraction U and for different values of the electron-phonon 
coupling constant. The magnitude of the parameters has 
been chosen as to allow a comparison with existing NRG 
calculations . V = 0.0f59, ujo = 0.05, U = 0.1 and 
A = 0, 0.03, 0.07, 0.09 (corresponding respectively to continu- 
ous, dashed, dotted and dash-dotted lines). The level position 
corresponds to an electron-hole symmetric case. 



Hamiltonian of Eq. pTjl . The parameters values have 
been chosen to allow comparison with NRG results for an 
electron-hole symmetric casei^i^^ and illustrate the tran- 
sition from the Kondo to the strong polaronic regime. 
For A = we find the typical spectrum with a Kondo 
resonance at the Fermi level, in good agreement with^^. 
When increasing A, the effective Coulomb repulsion U 
decreases and the Fermi-level resonance slightly broad- 
ens while the first phonon sub-bands appear (see curve 
for A — 0.03 and U — 0.064 in the figure). For larger 
values of A, U becomes negative, the Fermi level reso- 
nance eventually collapses, and the system experiments 
a transition to a polaronic regime with an energy gap in 
the spectral density. Strictly speaking, the Friedel sum 
rule is still satisfied, with a delta-like peak at the Fermi 
energy which at the scale of the figure cannot be resolved. 
These results are in good agreement with those of NRG 
calculationsi^iii^. 



We have presented an interpolative approach for cal- 
culating the electronic properties in a nanoscale sys- 
tem coupled to metallic electrodes and exhibiting both 
electron-electron and electron-phonon interactions tak- 
ing the Anderson-Holstein model as a test case. This 
approach provides a generalization of the interpolative 
self-energy method^^i^iiS^ which has been applied suc- 
cessfully to many different problems— iS£iS^i22ii2£. 

In a first step we have considered the spinless 
Anderson-Holstein model and discussed in detail for this 
case how the interpolative argument can be extended 
to deal with electron-phonon interactions. We have 
checked the accuracy of the approximation by applying 
the method to the case in which the electrode is repre- 
sented by a finite chain of three discrete levels which can 
be solved exactly by numerical diagonalization. In par- 
ticular we have found that the method gives a very good 
description of the spectral density at low energies. In the 
continuous case the method describes correctly the expo- 
nential decay of the width of the Fermi level resonance 
with increasing electron-phonon coupling. 

In a second step we have analyzed the full spin de- 
pendent Anderson-Holstein model and discussed how 
to generalize the interpolative self-energy for this case. 
This scheme naturally recovers the results of the purely 
electronic interpolation method for vanishing electron- 
phonon coupling. The accuracy of the method for this 
case is checked both for a discrete electrode model and by 
comparison with existing NRG results for the continuous 
modeliS,'^. 

We conclude that the interpolative technique is very 
reliable and provides a powerful tool to explore the role 
of strong electron-phonon interactions in more general 
situations of current experimental interest like electronic 
transport through single molecules. Work along these 
lines is under progress. 
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APPENDIX A 



The results presented so far in this section have con- 
centrated in the case of electron-hole symmetry. The 
interpolative ansatz is however valid as in the spinless 
case for a general asymmetric situation. In the extreme 
limit of a fully empty or fully occupied level, the model 
becomes in fact equivalent to a spinless case as it is ev- 
ident comparing Eq. with Eq. ^ for {n^) — 0, 
{{fiahg) = 0) or {ha) = 1, {(haha) = !)■ Therefore the 
ansatz tends also to the exact solution in that limit. 



For a purely electronic case the correlation function 
{fi^fii) can be calculated from the one-electron Green 
function of the localized level^^. A similar procedure for 
the Anderson-Holstein model yields the following equa- 
tion: 

U < hcrha > +A < h„{P + b) >= -^^j'duj Y.a{uj)Gaiuj) 

(Al) 
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where the propagators are assumed to be causal and the 
contour integration is on the upper complex plane. 

This equation could be used to calculate {hafis) if the 
averaged value {na{b^ + b)) can be obtained indepen- 
dently. This can be provided by the equation of motion 
of the phonon propagator. 

a 

where the Green function appearing on the right hand 
side can be used to calculate the average {n„{b^ + b)). 
From the equation of motion of this last Green function 
we havei^: 

<< ^na;b^+b »uj= \D'-°\uj) « ^n^;^ha- >>, 

(T a a 

The average value '^^naiw + b)) can be calculated 
by integrating the above equation: 

<5]ri.(St + 6)>=A^^i?(«)MXQH (A3) 
where xq(w) is the charge susceptibility appearing in Eq. 

El 

An interpolative argument can also be used to calcu- 
late xq{^^- First, we consider the atomic limit (F 0) 
and use the canonical transformation that diagonalizes 
the Anderson-Holstein Hamiltonian in this case. We ob- 
tain: 

X J*^ (co) = -27rz ^ [< > + < Kn-„ >] (A4) 

On the other hand, in the limit A, C/ ^ we can eval- 
uate this charge susceptibility function using Wick's the- 



orem: 

xg^-'=)(u;) = -27rz^ [< >2 + < >< h-, >] ^(c.) 

(7 

+ E^^°'('^) (A5) 

a 

where (uf) is the lowest order polarization bubble. In 

the hmit F ^ 0, ^a\uj) —2'Ki^h„){\ — {h„))b(uj) and 
Eq. lA5l tends to: 

x[]^''''\w) = -27^^E[< h„> ^ <h„><hs >] 

(A6) 

Due to these properties, the function: 
Xq{^) = ~27riE [< >^ -)- < fiafis >] (5(tj)+y^ Xo-(^) 

(J a 

(A7) 

interpolates correctly between both limits, Xcr{^) being 
the full polarization bubble. We have checked that a 
simple effective mass approximation in which Xcr(w) is 
calculated by 

X.(c.) = j dTGlff{T)Glff{~T)e-^-^ (A8) 

where C^J^ is a renormalized propagator with V^^^ = 
V/ (1 — dY.{^)/dujy^^, gives a good estimation of {fiafi^) 
for a broad range of parameters. As an additional check 
{fiahg) has been obtained by means of clusters calcula- 
tions with a good agreement with the above approxima- 
tion. 
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